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We develop a conceptually simple scheme based on a master-equation approach to evaluate the 
full-counting statistics (FCS) of elastic and inelastic off-resonant tunneling (cotunneling) in quantum 
dots (QDs) and molecules. We demonstrate the method by showing that it reproduces known results 
for the FCS and shot noise in the cotunneling regime. For a QD with an excited state, we obtain 
an analytic expression for the cumulant generating function (CGF) taking into account elastic and 
inelastic cotunneling. From the CGF we find that the shot noise above the inelastic threshold in the 
cotunneling regime is inherently super-Poissonian when external relaxation is weak. Furthermore, a 
complete picture of the shot noise across the different transport regimes is given. In the case where 
the excited state is a blocking state, strongly enhanced shot noise is predicted both in the resonant 
and cotunneling regimes. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.-b, 05.40.-a 


I. INTRODUCTION 

The full counting statisticJi^ (FCS) of charge transfer in 
quantum dots (QDs), nanostructures and molecules is an 
important component in the characterization of the mi¬ 
croscopic processes governing the transport. As the FCS 
contains the full information about the low-frequency 
current fluctuations, it provides access to all higher-order 
moments of the current fluctuations and hence gives in¬ 
sight not encoded in the average value of the current and 
the shot nois^ given by the two first moments. 

Experimentally, measurements of the FCS have been 
realized via real-time detection of single-electron tun¬ 
neling events,and the FCS of the charge-transfer 
mechanisms in various conductors ranging from tun¬ 
nel junctions and resonant tunneling in Coulomb block¬ 
aded QDs to Andreev tunneling at a superconductor- 
normal-metal interface has been characterized.^^^^ The¬ 
oretical sche mes to evaluate the FCS in different trans¬ 
port regimes pUMUl have successfully explained the ex¬ 
perimentally measured FCS as well as the observa¬ 
tion of super-Poissonian shot noise in resonant tunnel¬ 
ing through QDs.li^ Recent theoretical and experimen¬ 
tal work has focused on non-Markovia n effe cts due to 
coupling to external equilibrium bath^l^ and quan¬ 
tum coherent effects,^ finite-frequency current statis- 
ticsjinili] the effect of the electron-phonon interaction 
on the FCS in molecular contacts,l22H2ll interaction ef- 
fects,l22H2ll a,s well as the signature of Majorana bound 
states in the FCS.^^ 

Transport in the off-resonant regime where the QD 
levels are located outside the bias window and sepa¬ 
rated from the chemical potentials of the leads by an 
energy 5, is dominated by cotunneling processes.^^In co¬ 
tunneling processes, an electron or hole tunnels, either 
elastically or inelastically, through the energetically for¬ 
bidden charge state and only occupies the QD virtually. 


Inelastic cotunneling processes turn on at bias voltages 
exceeding the energy A of excited QD states and leave 
a clear fingerprint in the current-voltage characteristics. 
Such cotunneling spectroscopy is ideal for probing ex- 
cited states and their lifetime in QDs,^ molecules 
and grapheneQDs.l^SHUl Recently, the study of energy 
dissipatioiP^Illl and heat transporlP^ in the cotunneling 
regime has gained interest. Experimentally, the shot 
noise in the cotunnel ing re gime has been demonstrated to 
be super-Poissonian, HSHSi] in agreement with theoretical 
predictions.!^ Other theoretical studies have addressed 
the shot noise in the presence of cotunneling in specihc 
system^sHHHand the sig nature of cotunneling-assisted se¬ 
quential tunnelin j^^F7 | jn the shot noise.!^ 

The evaluation of the FCS taking into account co¬ 
tunneling, has been addressed theoretically.!^^!^ How¬ 
ever, the FCS in the cotunneling regime where the in¬ 
terplay between elastic and inelastic cotunneling governs 
the FCS, remains unexplored. Due to the difficulty of 
measuring single cotunneling events without collapsing 
the virtual intermediate state of a cotunneling process,!^ 
FCS of cotunneling processes has not been studied ex¬ 
perimentally. Recent theoretic al pr oposals for probing 
the cotunneling time Tcot ~ ;^/^j 63 | 64 ]may prove useful in 
that regard. 

In this paper, we develop a conceptually simple scheme 
based on a Markovian master equation description to 
evaluate the FCS of cotunneling processes. Compared to 
rigorous perturbative evaluations of shot noise and FCS 
to next-to-l eading o rder in the tunnel coupling strength 
F = 27rp|tp,!^2llSEn! the approach outlined here does not 
account for renormalization and broadening of the elec¬ 
tronic levels due to quantum fluctuations which give rise 
to non-Markovian dynamics.!^ However, for eV ^ 
F and in the cotunneling regime <5 3> F, non-Markovian 
effects are suppressed and can be safely neglected.!^ Our 
approach applies in these regimes, and we show that it 
recovers results for the shot noise and FCS in the cotun- 
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neling regime of simple models obtained with methods 
taking into account non-Markovian effectsP^^^ 

We furthermore demonstrate the method by studying 
the shot noise across transport regimes in a generic model 
for a QD system with an excited electronic state. In 
particular, we address the signature in the shot noise 
of cotunneling-related transport channels as well as the 
impact on the shot noise in the case where the excited 
state is a so-called blocking state. 


II. MASTER EQUATION APPROACH TO 
QUANTUM TRANSPORT 

Quantum transport in QDs and molecules involving 
higher-order tunneling processes between the QD and the 
leads, can be described with a T-matrix based master 
equation approach.l^In this approach, the rates for tun¬ 
neling between the QD states are evaluated using a gen¬ 
eralized Fermi’s golden rule by expanding the T matrix to 
a desired order in tunnel-coupling Hamiltonian, Ht- To 
lowest and next-to-leading order, this gives rise to sequen¬ 
tial and cotunneling processes, respectively. Compared 
to rigorous density-matrix approaches, where a formally 
exact master equation for the reduced density matrix of 
the QD can be obtained and systematically expanded 
in Ht (see, e.g., Refs. the T-matrix approach 

does not account for quantum effects such as broaden¬ 
ing and renormalization of the QD levels. It is therefore 
only applicable in the regime fc^T, eU ^ F as well as in 
the Coulomb blockade regime (cotunneling regime), i.e. 
i5 ^ F. In these regimes, the discrepancy between the T- 
matrix approach with proper regularized rates and exact 
perturbation theory vanishes.l^ 

In the T-matrix approach, the master equation is re¬ 
stricted to the diagonal elements of the reduced density 
matrix which are equivalent to the occupation probabil¬ 
ities Pm for the QD states—here labeled by a collective 
index m = (IV, i) for charge and excited state (electronic, 
vibrational, ...). The master equation governing their 
time evolution is given by 

Pm ~ Pm ^ ( ^m,m' T ^ ( Pm'^m^m- (1) 
m'^m m'^m 


Together with the normalization condition J2mPm = Ij 
it can be solved for the steady-state occupation proba¬ 
bilities, Pm = 0. 

Without the normalization condition, the master equa¬ 
tion takes the form of the matrix equation 


The transition rates Tmm' due to tunneling are ob¬ 
tained from the generalized Fermi golden rule 

9_ 

Tmm' = J-Y.\{nm?5{Ef-E,). (3) 

i'f> 


Here, \i/f) = \mlm') 0 \i'If) are products of QD and 
lead states, the sum is over possible initial \i') and hnal 
I/') states of the leads, and T = Ht-\-HtGoHt-\-- ■ ■ is the 
T-matrix with Gq = e -Ho denoting the Green function 
of the decoupled QD and leads described by Hq = Hqt, 
a = L,R, Ho, = and Ht = 

'Yjoika h.C.. 

The lowest-order sequential tunneling processes con¬ 
nect neighboring charge states of the QD system. In this 
case, the master equation takes the form 


PN,i Iggq 



+ r 





Pn-\-i,T°‘ N,i 

N + l,j 


+ PAf -1 jT“iv,i 
N-i,i 


(4) 


In next-to-leading order, cotunneling processes involve 
tunneling in and out of two leads and may excite the QD 
but do not change the charge state. These processes give 
rise to the following additional terms 

PN,i\cot = -PN,^ E +'^PN,3^f,p- (5) 


Note that the terms with z = j in the two sums cancel. 
These terms are associated with elastic cotunneling pro¬ 
cesses which do not change the state and therefore do not 
appear explicitly in the master equation. 

In addition to transport-induced transitions, relax¬ 
ation mechanisms due to coupling to bosonic degrees of 
freedom of an equilibrium environment (e.g., phonons) 
give rise to additional transitions between QD states. 
The transition rate for these processes is given by 

rZl> = '^\nB{AEmm')\, (6) 

where 7mm' determines the relaxation rate, AEmm' = 
Em' — Em, and is the Bose-Einstein distribution. In 
the absence of tunneling-induced transition, this results 
in a thermalized distribution of the QD states. 

From the steady-state solution, the current into ter¬ 
minal a can be obtained by evaluating the net rate of 
electrons. 


P = Mp, (2) 

where the diagonal (off-diagonal) elements of M (an 
M X M matrix where M is the total number of QD 
states) are given by the first (second) sum in Eq. Q. The 
master-equation matrix M is singular with the eigenvec¬ 
tor of the zero eigenvalue corresponding to the steady- 
state solution. 



where the first (second) term is the sequential (cotunnel¬ 
ing) current. 
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III. FULL COUNTING STATISTICS 

The main object of interest in counting statistics of 
charge transfer in QD systems is the probability distri¬ 
bution P{n, t) for n electrons having passed through the 
system from the source to the drain contact in the time 
interval t. 

In practice, it is more convenient to work with the 
cumulant generating function (CGF) S{x,t) which is de¬ 
fined by 


gS(x.t) (8) 

n 

where y is a counting field and S{0,t) = 0 in order to 
fulfill the normalization condition P{n,t) = 1. From 
the CGF, the cumulants of the current can be obtained 
as the derivatives with respect to the counting held y, 
i.e. lx->o for the m’th cumulant. The 

average current and the current noise are given by the 
hrst two cumulants. The probability distribution P{n,t) 
can be obtained by inverting Eq. Q, 

1 

P{n,t) = —j^ (9) 

which follows from the fact that S (x, t) is periodic in x 
with a period of 2tt. 

In order to calculate the CGF and the full counting 
statistics, it is convenient to work with the n-resolved 
probabilities, Prn{n,t), for the occupation of the states. 
As above, n here refers to the number of electrons hav¬ 
ing traversed the junction, and the distribution for the 
charge transfer is related to the n-resolved probabilities 
as P{n,t) = 

At the level of the master-equation treatment outlined 
in the preceding section, the time evolution of the n- 
resolved probabilities Pmin^t) is governed by 

=^M.{n-n')ip{n',t), (10) 

n' 

where the matrix elements of M describe the effect of 
tunneling and relaxation on the occupations Pm {n, t ), and 
it has been assumed that the dynamics is independent 
on the absolute value of the counting variable n and only 
depends on the difference n — n'. For sequential and co¬ 
tunneling processes, the change in the counting variable 


is restricted to the values n — n' = 0,±1, implying that 
Pm{ji,t) is only connected to the neighboring probabili¬ 
ties pra'{n±l,t) and Pm'^m{n,t). 

After Fourier transforming Eq. (10) to x space, the 
counting-field dependent master equation takes the form 

p(x,t) = M(x)p(x,f) (11) 


where 


P(X. 0 = X! (12) 

n 

and similarly for M. 

In the Markovian approximation, the GGF in the sta¬ 
tionary limit t —>■ oo can be obtained from the eigenvalue 
Amin(x) of tho counting- field d ependent matrix M(x) 
with the smallest real part,^!^^ i.e. 

‘5(x,l) = iAmin(x), (13) 

where t isthe measurement time. The evaluation of the 
FGS thus boils down to constructing the matrix M(x) 
for the Markovian master equation 0 and calculating 
the eigen value Aiiiin(x) from which the cumulants of the 
current can be obtained. 

The method developed by Bagrets and Nazarox®! ap¬ 
plies to sequential tunneling. In this case, the counting- 
field dependent master equation 0 can be constructed 
by replacing the rates in the second line of Eq. Q, which 
reside in the off diagonal of M, by the counting-field de¬ 
pendent rates,^ 


r“ ,fYj=r“ 1141 

where ± is for processes into/out of the counting lead a. 

Below we generalize this approach to cotunneling by 
demonstrating how to construct the X"dependent matrix 
M(x) when cotunneling processes are included. Our ap¬ 
proach is valid in the regime where ksT, eV ^ F or 
(5 3> F. In other words, when the tunneling-induced 
broadening F is smaller than one of the other energy 
scales. 


A. X'dependent master equation for cotunneling 

In order to derive the X"dependent master equation 
taking into account cotunneling processes, we start by 
writing up the master equation for the n-resolved prob¬ 
abilities, 


Pm(n,t)|^„t =-Pm(n,t) ^ F 


Q:/3 

mm' 


1 - SasS, 


OLO^mm' 


V 

/ . ^ m'r. 


(1 ^mm' ^Pm' iP't “1“ ^olL^^RP m' 1, t) + ^ctR^^LPm' “1“ 


(15) 
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where the terms with m = m! (m ^ m!) correspond to 
elastic (inelastic) cotunneling processes and we are count¬ 
ing the number of electrons n collected in the left lead. 


By differentiating Eq. (121 with respect to time and 
using the identity, 


^ e*^>(n ±l,t) = t), 


(16) 


we find 

Pm{X,t)\cot = 

m' ,oc^ 


1 


T 

/ ^ rn'r 




Pm'{X,t) 


= PmiX,t) 

m' ,cx0 


5mm' {5al3 + ^aL^pRC^^ + SaRSpLe — 1 


-f y 

' / ^ m'r 

= Pm{X,t) 

m' 


5a/3 + SaL^ISRe^^ + SaRSpLC 


Pm'iX,t) 


s /r“^ Cvi — , 

^mm ^ rnmyX) ^ mm' 




+ ^m'mix)Pm'{X,t), 


m' ^m^Oi/S 


where the y-dependent rates are defined by 

(x) = { 6^0 + S^rS^rc^^ + S^RSf^RC -^^). 


(17) 


(18) 


Here, the Kronecker delta Sap in the first term ensures 
that the contribution from elastic cotunneling to the rate 
in the second term inside the square brackets of Eq. 0 
is canceled when a = /?. These terms correspond to 
elastic cotunneling processes involving only one lead and 
therefore do not affect the counting statistics for the cur¬ 
rent. 

The x-dependent master equation in Eq. defines 
the cotunneling contribution to the matrix elements of 
M(x) in Eq. (Ill, and is our main formal result. It 


demonstrates how cotunneling processes can be included 
on equal footing with sequential tunneling processes by 
introducing counting-field dependent cotunneling rates 
[Eq. (@] in the master equation. Contrary to the con¬ 
ventional master equation for cotunneling [Eq. ([^], the 
X-dependent equation 0 contains contributions from 
elastic cotunneling processes [first term in the last line 
of Eq. 0]. This is due to the fact that, while they do 
not change the state of QD, they contribute to the trans¬ 
port and hence affect the charge-transfer statistics. We 
end by noting that the X'dependent elastic cotunneling 
rates appear in the diagonal elements of M(x), while the 
X-dependent inelastic cotunneling rates appear in the off 
diagonal. 


B. FCS of cotunneling 

In the following subsections, we show that our ap- 
proa ch re covers known results for the and shot 

nois^miSl pf cotunneling in simple models. 


1. Elastic cotunneling through a single level 

We start by considering the simple example of elastic 
cotunneling through a single off-resonant electronic level 
with energy eg- In this case, the X“dependent master 
equation 0 reduces to 


Poix,t) =Poix,t) 


where 


_ 

^ 00 ~ 


BoT^ 

2nh 


oijS 


de ■ 


ro“o^(x)-ri 


= IW(x)po(x,t), 

(19) 


, ^/o(e)[l-//3(e)] (20) 

(e - eo) 

is the rate for elastic cotunneling through the level and 
Bo(f (x) is the corresponding x-dependent rate given in 
Eq. (pj). 


The x-dependent matrix M(x) is here a scalar, and 
the eigenvalue Aniin(x)i nnd hence the CGF, can be read 
off directly from Eq. 


5(x) = to - 1) + ro“(e-*>^ - 1)] . (21) 

This CGF corresponds to bidirectional Poisson statistics, 
and is in agreement with previous work on FCS of elastic 
cotunneling through a single level [see Eq. (7) of Ref. 155] . 
For the current and noise we find 


I = 


dS 


to d{ix) 
BaT^ / 

2Trh I 


_ -pLR -pRL 

~ 00 ^00 


X^O 

de ■ 


1 


(e - eo)' 


IfUe) - /i^(£)] 


( 22 ) 
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and 


S = 


d'^S 


to diixY 

= coth 


= r, 


LR 

00 


X=0 


/ ey 


2Trh 


^RL 

00 


de 


{s-SoY 


[hie) - /fl(e)], 

(23) 


respectively, where eV = hl — t^R- At finite temper¬ 
ature, the Fano factor F = S/e\I\ is given by F = 
coth(eF/2A:BT). In the limit ksT » el^, the noise 
reduces to the equilibrium Johnson-Nyquist noise S = 
2GdkBT where Gd is the conductance. For eV ksT, 
shot noise becomes dominant with a Poissonian Fano fac¬ 
tor F = 1 as expected for independent tunneling pro¬ 
cesses. 


2. Elastic and inelastic cotunneling through a two-level 
system 


where Fy = F-®^ + is the total transition rate 

due to relaxation and inelastic cotunneling. 

To obtain the CGF, we set up the counting-field 
dependent master-equation matrix following Eqs. 0 
and (18), 


M(x) = 


roo(x)-roi rio(x) 
roi(x) rii(x)-rio 


(25) 


where the counting-field dependent rates are defined as 

r«(x) = rf«- 1 ) + - 1 ), ( 26 ) 


and 


r., (X) = rjf + Ff/ + + e*^F“ + (27) 

respectively. Note that rii(x = 0) = 0, imp lying that the 
standard master-equation matrix [Eq. & is recovered 
for X = 0. The derivatives of the counting-field depen¬ 
dent rates, which will be needed below, are given by 


Next, we consider cotunneling trough a system which, 
in addition to its ground state |0), has an excited elec¬ 
tronic state |1) with energy A relative to the ground 
state. In this case, inelastic cotunneling processes in¬ 
duce transitions between the ground and excited states. 
The steady-state occupation probabilities of the states 
are given by (see App. 0 


Po = 


10 


Gio + r. 


and Pi = 


01 


01 


10 


01 


(24) 


gr(x) 

d{ix) 

d^Tjx) 

d{ixY 


= Tlr^^ - 


(28) 

(29) 


for both the elastic and inelastic cotunneling rates. 

From the relevant eigenvalue of the counting-field de¬ 
pendent master-equation matrix (25), the CGF is found 
to be 


.(O ^ I 


roo(x) + rii(x) — Tio — For -f — Fio — roo(x) + rii(x)]^ + 4Foi(x)rio(x) 


(30) 


This CFG is a new result and describes the FCS of com¬ 
bined elastic and inelastic cotunneling. The FCS can be 
interpreted as arising due to switching between different 
bidirectional Poisson statistics, as discussed further be¬ 
low. 

The current and shot noise can be obtained from the y- 
derivatives of the CGF. In agreement with the standard 
master-equation calculation in App. we find for the 
current 


e dS 
io d{ix) 


x=o 


— Po^ix^oo + Pldix^ll + Podix^Ql + Pl^i^^-Fio 
= Ie\ + Anel, (31) 


where /ei and /inei are elastic and inelastic contributions 
given by the two first and two last terms in the second 
line, respectively. 

The noise given by the second derivative of the CGF 


is found to be 


S = 


d^s 


to d{ixy 


x=o 


— Po{dix^oo + dj^oi) +Pi{d^x^ii + d'^iXio) 

+ -—p— 5ip(.roii9i;^Fio — ^ix^oo^ixrii 

r 01 + r 10 L 

+ (Sij^Foo + iSixTii) I — 

= 5'poisson + AF, (32) 


where the two terms in the second line describe equi¬ 
librium Johnson-Nyquist and Poissonian shot noise and 
the term in the square brackets is responsible for a non- 
Poissonian correction AS” at bias voltages larger than 
the inelastic threshold E > A. We note that this re¬ 
sult for the shot noise is in agreement with Ref. [5H [their 
Eq. (5.5)]. 
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Due to the factor in front of the square brackets, the 
non-Poissonian correction AS” diverges for PoijPio —t 0 
if, at the same time, Poo ^ Pii- The diverging super- 
Poissonian noise can be understood as follows. In the 
limit roi,rio — t 0, i.e. for negligible environmental re¬ 
laxation and vanishing inelastic cotunneling rates, inelas¬ 
tic cotunneling processes change the state of the system 
at a rate that is slow compared to the rate of elastic co¬ 
tunneling processes. As the latter dominates the current, 
this results in a current that switches between different 
values when Tgo ^ Pn- Such telegraphic switching be¬ 
tween two transport channels with different conductance 
naturally produces super-Poissonian shot noise. 

3. Quantum dot with a spin-split level 

The shot noise in a QD with a spin-split level has 
previously been studied with the real-time diagrammatic 
method in Ref. This goes beyond the T-matrix ap¬ 
proach adopted here by accounting for the broadening 
and renormalization of the QD levels. Here, we show 
that the two methods yield identical results for the noise 
in the regime ksT ^ P. For ksT ^ P, the disagreement 
between the two methods amounts to a small quantita¬ 
tive difference. 

The Hamiltonian of the QD is given by 

Hqu = -b Un^ni, (33) 

cr 

where Sa- is the spin-dependent level position and U is the 
Coulomb interaction for double occupancy of the QD. 
The sequential and cotunneling rates are calculated as 
outlined in detail in Sec. IIV Al below. 

In Fig.[^we show the Fano factor obtained with our ap¬ 
proach as a function of bias voltage and for different val¬ 
ues of the ratio T/ksT at hxed temperature, ksT = 0.1. 
With the parameters specified in the caption of Fig. 
the spin-down level is filled at H = 0. The inelastic spin- 
flip channel opens at eV = A = — ei- At eH/2 = 

(eH/2 = U), the removal (addition) energy for 

the spin-down (spin-up) level becomes resonant with the 
chemical potentials of the drain (source) electrode. 

As mentioned above, the agreement with the real¬ 
time diagrammatic approach of Ref. is perfect for 
ksT ^ F. Remarkably the agreement is even quite good 
for ksT > P (for the present parameters less than 10% in 
the whole bias range for P/fc^T = 0.5). For small values 
of F, the current and noise are dominated by thermally 
activated sequential tunneling through the two spin levels 
giving rise to super-Poissonian noise.l^ At larger values 
of P, the change in the Fano factor is due to cotunneling 
processes. At low bias, elastic cotunneling dominates and 
results in F = 1. Above the inelastic cotunneling thresh¬ 
old at eV = A, the state with an electron in the spin-up 
level becomes populated. After a cotunneling-induced 
spin flip (i—^t), tbe fact that the spin-up level is located 
inside the bias window, < eV, opens for transport 



FIG. 1: Fano factor for a QD with a spin-split level for 
different values of the ratio T/ksT near the cross over be¬ 
tween the cotunneling and resonance regimes (see Fig. 3 
in Ref. 1531 for comparison). Parameters: = —1.5, = 

0.5, C/ = 4, keT = O.I, Tl/r = F/2, ptR/R = ±eV/2. 

through the QD via sequential tunneling. This trans¬ 
port channel remains open until the QD relaxes to the 
spin-down state via a sequential tunneling or an inelas¬ 
tic cotunneling processes. The transport thus switches 
between being dominated by elastic cotunneling through 
the I level and sequential tunneling through the f level. 
As a consequence of this alternating change in the trans¬ 
port mechanism, the noise becomes super Poissonian. 


IV. QUANTUM DOTS AND MOLECULES 
WITH EXCITED ELECTRONIC STATES 

In this section, we consider the shot noise of a QD 
with an excited electronic state. This could, for exam¬ 
ple, be a low-energy QD level, a molecular orbital or spin 
excitations (spin manifolds). We here give a complete 
picture of the shot noise across the different transport 
regimes and study, e.g., the signature of the interplay 
between cotunnelin g and cotunneling-assisted sequential 
tunneling (COSETp^IIIl in the shot noise. We further¬ 
more consider the situation where the excited state is a 
so-called blocking state. This situation is familiar from, 
e.g., quantum dots and molecules with excited states or 
broken degeneracies.^sMH] jn the resonance regime outside 
the Coulomb blockaded regions, blocking states give rise 
to pronounced negative differential conductance (NDC) 
and strong super-Poissonian noise.^ In the cotunneling 
and COSET regimes, the effect of the blocking state on 
the shot noise has, so far, not been studied, and we find 
that the noise changes qualitatively in the presence of the 
blocking state. 


A. Generic model and transition rates 

We consider a spinless model for a QD where the N 
electron configuration has an excited electronic state. 
The states of the QD are described by a set of generic 
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many-body states 

I7V-1), |7Va), |iV6), |iV + l), (34) 

where N refers to the number of electrons on the QD. 
In a microscopic description of the QD, the states and 
their energies result from the diagonalization of the un¬ 
derlying microscopic Hamiltonian. The latter are here 
parametrized as 

En-1 = 0, Epfa{b) = So (+A) 
and En+i = 2eo + U. (35) 

where A is the energy of the excited A-particle state 
|iV6) relative to the ground state \Na) with energy fo = 
So — eVg (relative to the N — 1 state), and U is Coulomb 
energy associated with the addition of an electron to the 
fV-particle state. With this parametrization, the addi¬ 
tion and removal energies of the A-particle ground state 
become E^+i — Ej^a = So + U and Ej^a — = Soi 

respectively, implying that the transport gap of the QD 
is given by U. 

To take into account sequential and cotunneling pro¬ 
cesses, we expand the T matrix to second order in Hx in 
the calculation of the transition rates. 

The sequential tunneling rates for adding and remov¬ 


ing an electron from the QD are given by 

^N-l,Na/b = ~^\^N-l,Na/b\'^ fa{Sa/b) (36) 

^Na/b,N+l = ~^\^Na/b,N+l\‘^faiSa/b + U) (37) 

and 

^Na/b,N-l = ~^\^N-l,Na/b\'^ [l ~ fa(,£a/b)] (38) 

r^+l.ATa/b = ^ma/b,N+l\" [l ” Us^/b + U)] , 

(39) 

respectively, where Tq, = 2TTpa\ta\'^ is the tunnel broad¬ 
ening, ea(b) = So (+^)i and the matrix elements between 
the many-body states are given by 

^Ni,N-i = {N — l|d^|iVz) = (40) 

= {N+ l\dl\m) = (41) 

Here, dl^,d^ denote the creation and annihilation opera¬ 
tors for the single-particle states in the QD system that 
couple to lead a. For a given microscopic model, the 
matrix elements can be obtained from the many-body 
states. Here we treat them as tuneable parameters. 

The elastic cotunneling rates for the different states 
are given by 


r«/3 

FaT/J 

1 de 


^ N-1 

2Trh j 

CH 

1 

o 


FaF/S 

j de 

^^^Na,N-l^^^N-l,Na 

^ N,a 

2TTh J 

£ - £0 

-pa^ 

TaT/S 

f de 

-^^^Nb,N-l^^^N-l,Nb 

^ N,b — 

27rS j 

< 

1 

■ o 

1 

pa/3 

FqF^ 

f de 

M°‘ 

^'^^N+l,Na^''^Na,N+l 

^ iV+1 

2T:h j 

e — io — U 


-b 


^''^N-l,Nb^'^^Nb,N-l 

e-So- A 

^^^Na,N+E'^^N+l,Na 

e — e{) — U 

N-\-l,Nb 

E — S'o — tl ^ 

^^^N-El,Nb^^-^Nb,N+l 

£ — S'q — JJ ^ 


fa{s) [1 - ffi{s)] 

(42) 

) 

/a(e) [1 - //3(e)] 

(43) 

Uis) [1 - //3(e)] 

(44) 

/a(e) [1 - //3(e)]. 

(45) 


The inelastic cotunneling rates between the ground and excited state of the A^-electron configuration are given by 


N,ab 


NM 


FaFft 

1 de 

^^^Nb,N-l-^^^N-l,Na 

^^^Nb,N-\-l^^^N+l,Na 

2Trh J 

e - £o - A 

e — io — U 

TaTb 

jde 



2Trh J 

e - £o 

£ — £q — U ^ 


u{s) [1 - fpis-A)] 


U{s) [1 - fpie + A)]. 


(46) 

(47) 


We evaluate the cotunneling rates at finite tempera- scheme described in App. 
ture and bias with the commonly applied regularization 

A situation which resembles the conditions for strong 



































FIG. 2: (Color online) Stability diagrams showing the differential conductance G = dl/dV (upper row) and the 
Fano factor F = 5'/e|/| with the thermal contribution to the noise subtracted (lower row) as a function of gate and 
source-drain voltage for different situations. Note the different color scales in the lower row. (a),(b) The ground and 
excited fV-particle states have identical matrix elements to the A^± 1 states, = 1. (c),(d) and (e),(f) The 

excited state is a blocking state with small matrix elements to the ± 1 states: (c),(d) = 0.1, and (e),(f) 

^NbN+i — O-l- Matrix elements not specified are equal to unity. Parameters (in units of U): £o = ~l/2, U = 1, 
A = b.l, Fi/^ = 0.001, keT = 0.01, = ±1^/2. 


super-Poissonian noise discussed below Eq. (32), is real¬ 
ized if the excited state \Nb) is a blocking state which 
is characterized by having small matrix elements to the 
other states, i.e. <C 1. This leads to strongly 

reduced cotunneling rates and implies that the inelastic 
rates and the elastic rate for \Nb) are reduced compared 
to the elastic cotunneling rate for state |Aa). Strong 
super-Poissonian noise is therefore expected. Since the 
matrix elements to both the TV ± 1 states have to be 
small in order to suppress the elastic cotunneling rate for 
\Nb), the shot noise is highly sensitive to the properties 
of the blocking state via its matrix elements with the 
neighboring charge states. 


B. Stability diagrams and shot noise 

We now study the current, conductance, shot noise, 
and Fano factor for different situations for the blocking 
property of the excited state \Nb). In the upper row of 
Fig. [2] we show the differential conductance as a function 
of gate and source-drain bias voltage—also referred to as 
charge-stability diagrams—for the cases without [^a)] 
and with j^c) and[^^ a blocking state. For the latter, 
the two plots in Figs.|^c) and[^e) correspond to differ¬ 
ent situations for the matrix elements involving \Nb) (see 
caption of Fig. for details). Inside the Coulomb block¬ 
aded regions where sequential tunneling is suppressed. 


the current is dominated by cotunneling processes. Due 
to the linear color scale in the figures, the cotunneling 
features in the conductance are not visible. 

Outside the blockaded regions where sequential tunnel¬ 
ing dominate the current, excitation lines going out from 
the central blockaded region appear at voltages where the 
excited state enters the bias window. Some of these lines 
show pronounced NDC when the excited state is a block¬ 
ing state. Depending on details of the matrix elements 
for the blocking state, NDC occurs either for both signs 
[^c)] or one sign [^e)] of the gate and bias voltages, and 
may completely suppress the current He)]. 

The corresponding stability diagrams for the Fano fac¬ 
tor are shown in the lower row of Fig. Like the con¬ 
ductance, the Fano factors are highly sensitive to the 
matrix elements of the blocking state. The Fano factors 
have pronounced features with strong super-Poissonian 
values in the blockaded regions of the stability diagrams 
where the signature in the current and conductance is 
weak. These features originate from the opening of co¬ 
tunneling related transport channels. In particular, at 
the threshold for inelastic cotunneling at D = A, the 
Fano factor increases markedly. Also, at the onset of 
COSET processes near the edges of the blockaded re¬ 
gion, the Fano factor shows drastic changes. In the re¬ 
gions outside the blockaded regions where NDC occurs, 
strong super-Poissonian noise with a large Fano factor is 
observed. The Fano factor in Fig.j^f) where the excited 
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(a) (c) (e) 








FIG. 3: (Color online) Super-Poissonian noise in the cotunneling regime, (a) Rates for inelastic cotunneling between 
the \Na), \Nb) states as a function of bias, (b) Occupation probabilities for the \Na/b) states. (c),(d) Cotunneling 
current [(c)], shot noise [(d) left axis] and Fano factor [(d) right axis] vs bias voltage. (e),(f) Same as in (c),(d) for the 
situation where the excited state is a blocking state. The plots correspond to the low-bias part of the cuts indicated 
by dashed lines in Fig. j^b) [(a)-(d)] and Fig. [^d) [(e)-l-(f)]. 


state is only partially blocked, shows a mixture of the fea¬ 
tures present in Figs, [^b) and[^d) for no and complete 
blocking, respectively. 

In the following subsection, we analyze the shot noise 
in the different transport regimes in closer detail. 


1. Super-Poissonian noise in the cotunneling regime 


We start by considering the noise in the cotunneling 
regime (low-bias part of the cuts through the center of 
the diamonds in Fig.[^. In this regime, the results based 
on a pure cotunneling description from Sec. Ill B 2| apply 
[Eqs. ( |M| ) and j^]. 

Figure [^a)-j3p) summarize the situation without a 
blocking state corresponding to the stability diagrams in 
Figs, [^a) andj^b). The inelastic cotunneling rates are 
shown in Fig. I^a) as a function of bias voltage for the 
different lead indices. At low temperature, inelastic co¬ 
tunneling processes in the direction of the voltage drop 
(L —>■ R) may excite the QD when V > A. At the 
threshold P = A, the corresponding rate F^^^ starts 
to increase linearly with the applied bias. De-excitation 
processes with rates which relax the QD back to 

its ground state, are always possible. The resulting oc¬ 
cupation probabilities given by Eq. (241 are shown in 
Fig. I^b) with the probability for the ground (excited) 
state decreasing (increasing) approximately linearly with 
the bias near the threshold E > A. 


The occupation of the excited state at E > A gives rise 
to a strong inelastic signal in the current [Fig.j^c)] corre¬ 
sponding to a positive step in the differential conductance 
dl/dV. At high bias voltage E A, the elastic and in¬ 
elastic contributions to the total current in Eq. (31) be¬ 
come equal. The shot noise shown in Fig. I^d) (left axis) 
together with the Fano factor (right axisjf also shows a 
clear signal at the inelastic threshold. The Fano factor 
is given by its equilibrium value F ~ coth(eE/2A:BT) at 
low bias E < fc^T, drops to the Poissonian value F = 1 
for ksT < E < A, and becomes super-Poissonian with 
F > 1 for E > A. The modest value of the super- 
Poissonian Fano factor {F ^ 1.2) stems from the fact 
that the elastic and inelastic cotunneling rates are of the 
same order of magnitude. 


Figures l^e) and[^f) show the current and shot noise 
in the case where Nb) is a blocking state [cut along 
dashed line in Fig. ^d)]. In this case, the elastic and 
inelastic cotunneling rates involving \Nb) are strongly re¬ 
duced. Therefore, the current in Fig. [^e) is completely 
dominated by the elastic component through JAa), and 
the corresponding differential conductance has a nega¬ 
tive step at the inelastic threshold where the badly con¬ 
ducting excited state becomes populated. The noise in 
Fig. if), on the other hand, increases non-linearly with 
the bias for E > A. As a consequence, the Fano factor 
becomes strongly super-Poissonian with E 3> 1. The 
mechanism behind the increased noise can be under¬ 
stood as telegraphic switching between two elastic cur- 
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(a) (b) (a) 




5 

- F 

-- coth(eV/2fcj3T) 
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(b) 



FIG. 4: (Color online) Fano factor in the cotunneling and 
COSET regime along the cuts marked with dashed lines 
in the stability diagrams of Fig. The Fano factor F^o 
obtained from the pure cotunneling expressions for the 
current and noise in Eqs. © and ( [ 3 ^ is also shown. 
Above the onset of COSET processes a pure cotunneling 
description becomes ill defined. The plots correspond to 
the cuts indicated by dashed lines in Fig. I^b) [(a)] and 
Fig.[|;d) [(b)]. 


rent channels with different intrinsic shot noise. 

We emphasize that the appearance of strong super- 
Poissonian noise in the cotunneling regime relies on the 
excited iV-electron state being blocked from both the N± 
1 states, i.e. all the matrix elements must be 

small. In addition, the relaxation rate due to coupling 
to an external equilibrium bath in Eq. (§ must be small 
compared to the inelastic cotunneling rates, ^ 

^’tfab/ba- case, the QD relaxes to its 

ground state on a time scale much faster than the time 
between cotunneling events. Hence, elastic cotunneling 
via the ground state will dominate the current and noise 
and the Eano factor becomes Poissonian F ~ 1 in the 
limit of strong external relaxation. 


EIG. 5: (Color online) Eano factor in the resonant regime 
along the cuts marked with dotted lines in the stability 
diagrams of Fig. The plots correspond to the cuts 
indicated by dotted lines in Fig. I^b) [(a)] and Fig. I^d) 

[(b)]. 


COSET regime, the change in the Eano factor is quali¬ 
tative in the two cases. Whereas it increases to a value 
of F ^ 2 without, it drops to F = 1 with a blocking 
state. The increase in the Eano factor in the former case, 
is due to COSET processes where a channel for sequen¬ 
tial tunneling inside the blockaded region opens when 
the excited state becomes populated via inelastic cotun- 
neling.l^ This gives rise to a current that is alternately 
governed by sequential and cotunneling every time an 
inelastic cotunneling process excites and de-excites the 
QD, respectively. However, in the presence of the block¬ 
ing state, sequential and cotunneling via the excited state 
are strongly suppressed. The current is therefore dom¬ 
inated by elastic cotunneling via the ground state and 
the shot noise becomes Poissonian with F = 1. This is 
markedly different from the situation without a blocking 
state.^ At bias voltages V > 2|£o| and V > 2|£o + U\, 
the main resonances enter the bias window and sequential 
tunneling becomes dominant. 


2. Noise in the COSET regime 


To get a clearer picture of the behavior of the noise in 
the part of the Coulomb blockaded region where cotun¬ 
neling and sequential tunneling coexist and COSET pro¬ 
cesses provide a relaxation channel for the excited state, 
we show in Eig. the bias dependence of the Fano fac¬ 
tor along the full cuts (positive bias only) marked with 
dashed lines in Fig.[^ For comparison, the dotted lines in 
Fig. [ 4 ] show the Fano factor F^o obtained from Eqs. (31) 


and p2[ ) taking into account cotunneling only. At low 
bias, the noise is dominated by cotunneling processes. 
As discussed above, the noise becomes super-Poissonian 
at the threshold for inelastic cotunneling at P = A and 
acquires a strongly super-Poissonian Eano factor F ~ 15 
in the presence of the blocking state. 

The onset of COSET processes takes place at the side¬ 
band resonances at P = 2|£o + P —A| and P = 2|eo-l-A|, 
where relaxation of the excited state \Nb) via sequential 
tunneling to the \N ± 1) states becomes possible. In the 


3. Sub-Poissonian noise, NDC and super-Poissonian 
telegraphic noise in the resonant regime 

Figure shows the bias dependence of the Fano fac¬ 
tor along the cuts in the resonant regime outside the 
Coulomb blockaded regions marked with dashed lines in 
Fig.i In the situation without a blocking state, the Fano 
factor is sub-Poissonian with F = 0.5 for ksT < V < 2A 
and a slightly larger value F ~ 0.55 for P > 2A, which 
is characteristic for sequen tial tunneling through a QD 
with excited states.^ ^^F3 | However, in the presence of 
the blocking state, the shot noise becomes strongly super- 
Poissonian with F ~ 23. The mechanism responsible for 
NDC and the strong enhancement of the noise is the 
same. When the blocking state enters the bias win¬ 
dow, the QD gets trapped in the blocking state due 
to the small transition rate to other states. This re¬ 
duces the current and results in a telegraphic noise with 
long quiet periods without charge transfer interrupted 
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by avalanches of transfer processes every time the QD 
escapes the blocking state. 


be applied to investigate, e.g., the interplay between in¬ 
elastic cotunneling and quantum interference in the PCS 
of molecular contacts.l^ 


V. SUMMARY 


In summary, we have demonstrated how the standard 
scheme to evaluate the PCS of charge transfer due to se¬ 
quential tunneling in Coulomb blockaded QD systemJi^ 
can be generalized to take into account cotunneling pro¬ 
cesses. In analogy with the procedure for sequential tun¬ 
neling, this is done by replacing the cotunneling rates in 
the Markovian master equation Q with counting-field 
dependent rates as described in Eqs. (171 and ( 181. This 
approach neglects non-Markovian effect^SHHI^"^ associ¬ 
ated with tunneling-induced level broadening and shifts, 
and, hence, only applies for ksT, eV ^ P or in the co¬ 
tunneling regime, <5 3> P. In the cotunneling regime, we 
have demonstrated that the results for sh ot noise and the 
PCS from more elaborate method^^^Ml are reproduced. 
In addition, we have obtained an analytic expression for 
the CGP [Eq. (30l] describing the charge-transfer statis¬ 
tics of elastic and inelastic cotunneling in a two-state QD 
system. 

Studying a generic model for a QD with an excited 
electronic state, we found that the shot noise in the co¬ 
tunneling regime is inherently super-Poissonian for volt¬ 
ages larger than the inelastic threshold V > A. A 
strongly enhanced noise level with Pano factor F ^ 1 re¬ 
sults when the excited state is a so-called blocking state. 
This is due to telegraphic switching between the two dif¬ 
ferently conducting channels for elastic cotunneling via 
the ground and excited state. In the presence of environ¬ 
mental relaxation, the super-Poissonian noise is reduced 
and becomes Poissonian with F = 1 once the relaxation 
dominates the inelastic cotunneling rates. In the COSET 
regime where cotunneling and sequential tunneling co¬ 
exist, we found that the noise is, respectively, super- 
Poissonian and Poissonian for an excited state without 
and with blocking properties. 

Our approach for evaluating the PCS can be general¬ 
ized to other higher-order tunneli ng pr ocesses in QD sys¬ 
tems, such as, e.g., pair tunnelinJ^^IZIl and charge recon¬ 
figuration processes in multi-dot systemsand may 
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Appendix A: Master equation for elastic and 
inelastic cotunneling in a two-state system 

In this case, the standard master equations takes the 
form of the 2x2 matrix 

M = (d’” f;” ), (Ai) 

\ t 01 —t 10/ 

where Tij = rates for the inelastic 

cotunneling-induced transitions between the states. 

The eigenvalues are found to be 

A = - [— (Pio + Poi) ± (hio + Poi)], (A2) 

with the eigenvector of the zero eigenvalue giving the 
steady-state solution. 


f Tip Pqi \ 

'vPoi + rio’roi + Pioy 


(A3) 


The current can be obtained from the steady-state solu¬ 
tion as 


I = po (roY - Polf) + Pi - Pfi^) 

+ Po(ro^i^-ro“)+Pi (rfo«-r“), 

= del + dineh (-^4) 

where the first and second line are the elastic /ei and 
inelastic /inei contributions to the current, respectively. 


Appendix B: Regularized cotunneling rates 

The integrals for the elastic and inelastic cotunneling rates can be evaluated using the standard regularization 
scheme. The procedure for regularizing the diverging integrands in the cotunneling rates of Eqs. (42)“(47) can be 
found in Ref. IHU 

The integrants in the expressions for the cotunneling rates can be written on the general form 


A 


e — El 


i0+ 
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e- 62 


i0+ 
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£ — El 


■ 10 + 
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£ — £2 + *0+ 


±2Re 
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£ — £1 + 10+ £ — £2 — *0+ 


(Bl) 


where we have added a infinitesimal broadening of the QD states (the regularizer) in the denominators. 
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The regularized rates can then be obtained analytically from the following two integrals, 


de 


f{e - E,)[l - f{e - E2)] _ nB{E2 - E,) 


(e - ei)(£ - £ 2 ) 


£1 - £2 


Re ['tl^{E+^) - ip{E^ 2 ) - + ^(£' 12 )] 


(B2) 


and 


/ 


de 


/(^ 


i;^)[l-/(£-i^2)] 

(£ - £1)^ 


nB{E2 — El) 

2TrkBT 


Im - V^'(£^n)] , 


(B3) 


where ip denotes the digamma function, ip' its derivative, Ef^ = ^ ± ~ Bose-Einstein 

distribution. 
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